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Abstract. Radiative transfer models of multi-transition observations are used to determine molecular abundances as functions 
of position in pre- and protostellar cores. The data require a "drop" abundance profile with radius, with high abundances in the 
outermost regions probed by low excitation 3 mm lines, and much lower abundances at intermediate zones probed by higher 
frequency lines. The results are illustrated by detailed analysis of CO and HCO* lines for a subset of objects. We propose a 
scenario in which the molecules are frozen out in a region of the envelope where the temperature is low enough (< 40 K) 
to prevent immediate desorption, but where the density is high enough (>10'*-10' cm"') that the freeze-out timescales are 
shorter than the lifetime of the core. The size of the freeze-out zone is thereby a record of the thermal and dynamical evolution 
of the cores. Fits to CO data for a sample of 16 objects indicate that the size of the freeze-out zone decreases significantly 
between Class and I objects, explaining the variations in, for example, CO abundances with envelope masses. However, 
the corresponding timescales are 10^*""^ years, with no significant difference between Class and I objects. These timescales 
suggest that the dense pre-stellar phase with heavy depletions lasts only a short time, of order 10^ yr, in agreement with recent 
chemical-dynamical models. 
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1. Introduction 

The environments of the youngest pre- and protostellar ob- 
jects are characterized by large amounts of cold gas and 
_ dust. The chemistry in these early stages is affected to a 
■ large degree bv freeze-o ut of molecules onto dust grains (e.g., 
iBergin & Lang eii Il997 ^ and is closely related to the thermal 
evolution of the cores. An important question is how the pre- 
' and protostellar stages are linked and what their respective 
timescales are. In the pre-stellar stages the thermal balance is 
dominated by the external radiation field, which, except in spe- 
cial cases, does not heat the mater ial to temperatures higher 
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than w 15 K (e.g., Evans et al., 2001). At such low temperatures 

most molecules gradually freeze out, with very short timescales 
in the innermost dense regions and increasing timescales to- 
ward the exterior where they may become longer than the age 
of the core (e.g.. Caselli et al.. 1999). 

In the protostellar stages, in contrast, the thermal balance 
is dominated by the heating from the central newly formed 
protostar, introducing a steep temperature gradient toward the 
core center. Radiative transfer modeling of the dust continuum 



Send offprint requests to: Jes K. J0rgensen 

* Present address: Harvard-Smithsonian Center for Astrophysics, 
60 Garden St. MS42, Cambridge, MA 02138, USA 
Correspondence to: jjorgensen@cfa.harvard.edu 



emission shows that the characteristic temperatu res can rise to 
a few hundred K in the innermost regions ('e.g.. lShirlev et all 
2002; J0rgensen et al., 2002). Still, significant depletions are 
observed also in these stages, e.g., for CO, indicating that a sub- 
stantial fraction of the envelope material remains at low tem- 
peratures (e.g., Blake et al., 1995; Ceccarelli et al., 2001) . In a 
large survev o f pre- and protostellar objects, |j0rgensen et alJ 
1I2OO21 1200 41^ found a strong correlation between the abun- 
dances of CO (and related species such as HCO^) and envelope 
mass. 

The timescales for freeze-out and evaporation depend sen- 
sitively on density and temperature. The availability of accu- 
rate physical structures from dust continuum data thus provides 
an opportunity to constrain the timescales independently us- 
ing only chemistry. Currently, the ages of pre- and protostel- 
lar objects are determined almost exclusively from statistics, 
such as the number counts of cores wi th and withoti t asso- 
ciated far-infrared (IRAS) sources (e.g.. iLee"^ Mversill999l: 
J essop & Ward-Thompso n. 2000). resulting in a large spread in 
pre-stellar ages from ~ 10^ to a fewx 10'' yr. Our semi-empirical 
procedure for constraining the chemistry and its assumptions 
are summarized in Fig. 1 of .Dotv et al. (2004). and consists in 
the simple st case of fitting a con s tant ab undance to the data, 
as used in J0rgense n et al.l 1I2OO2', 2004b). In those papers, it 
was also realized that constant abundance models provide poor 
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fits to the lowest excitation lines of species such as CO and 
HCO^, leading to the proposal of a "drop" abundance profile 
with a specific freeze-out region over part of the envelope. Such 
"drop" profiles were also found to best reprod uce the inter fer- 
ometer data of H2CO in two sources ( Schoier et all l2004l) . In 
this paper, we further explore and quantify the "drop" abun- 
dance profiles for the full set of 16 objects and relate the size 
and location of the depletion zones to pre- and protostellar evo- 
lution. This semi-empirical study forms an important comple- 
ment to full chemo-hydrodynamical models which follow the 
che mistry in time as the matter collapses to form a central star 
(e.g. jRawlings et allll99ilLee et alll2004 . 

2. Model 

The main assumption in our analysis is that the chemical struc- 
ture in the protostellar stages is controlled by thermal desorp- 
tion processes. Although other non-thermal processes such as 
cosmic-ray induced desorption play a role for weakly-bound 
species like CO, they cannot prevent fre eze-out in the densest 
and coldest gas re.g.. lShenetani2004 . Similarly, photodes- 
orption can keep molecules off the grains, but this only involves 
the outermost regions up to an Ay «; 2 (Bergin et al., 1995). 
These effects are expected to be secondary to the abundance 
structure resulting from the freeze-out and thermal evaporation 
of CO. Finally, outflows and shocks may be important in regu- 
lating the abundance structures, but the narrow line-widths for 
optically thin species such as C'**0 and H'^CO^ suggest that 
they probe predominantly the quiescent bulk envelope mate- 
rial. Our analysis focusses on the chemistry in the outer enve- 
lope on 500-10,000 AU scales, and does not consider the addi- 
tional abundance jumps in the innermost region (<100 AU) at 

T > 100 K w here all ices evaporate . 

Following iRodgers & Charnlevl l Eool . the thermal desorp- 
tion rate ^ and the freeze-out rate A can be written as: 

^(M) oc exp(-^) (1) 
/ T \°-^ 

i(M) =4.55xlO-'M— f- «H [s-'] (2) 
\ m(M) / 

where and Tg are the dust and gas temperatures, respec- 
tively, m(M) the molecular weight, and hh the total hydrogen 
density. £'b(M) is the binding energy of the molecule depending 
on the ice mantle composition, for which we adopt the values 
tabulated by Aikawa etal,(1997J (£b=960 K for CO). The CO 
desorption temp eratures of 30 K or high er inferred from the ob- 
servational data ('j0rgensen et al.', 2002) suggest that not all CO 
is bound in a pure CO ice matrix but that at least some of it is in 
a mixture with H7O where the b inding energy may increase to 
-1200 K Pollin gs etalll2003l) . The pre-exponential factor in 
Eq. Q depends on whether "zeroth" or "first" order desorption 
kinetics are considered. In astrochemical models first-order ki- 
netics are commonly used, which is appropriate for the desorp- 
tion of (sub-)monolayer quantities of an adsorbate from a solid 
surface. For a thick layer of pure CO, a zeroeth-order formu- 
lation of the desorption process is more appropriate, however 
(see discussion in <Collings et al.„..2003.) . Our calculations adopt 
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Fig. 1. Comparison between CO desorption and freeze-out 
timescales as functions of temperature and density. The solid 
line indicates the desorption timescale while the dotted and 
dashed lines indicate the (density and temperature depen- 
dent) freeze-out timescales for TMRl and N1333-I2, respec- 
tively. Depletion occurs where the curves for the freeze-out 
timescale intersect the dark colored region, in this example 
for an assumed age of 10^ years (dash-dotted line). Freeze- 
out timescales corresponding to H2 densities of 1x10^, 1x10^ 
and 1x10^ cm^^ for NGC 1333-IRAS2 and 1x10^, 1x10^ and 
1x10^ cm"^ for TMRl have been indicated by the filled and 
open circles, respectively. 



a first-order formulation with a pre-exponential factor of 10'^ 
but this is not expected to change our main conclusions since 
the exponent stays the same and dominates the temperature be- 
havior. 

Fig. [2 shows the desorption and freeze-out timescales of 
CO, defined as 1/^(M) and 1/A(M), respectively, as functions 
of depth for two protostars from the sample of .Icir gensen et alJ 
(l2002h . NGC 1333-IRAS2 and TMRl . These objects are classi- 
fied as Class and I, and have significantly different envelope 
masses of 1.7 and 0.12 Mq, respectively. Their temperatures 
and densities have been constrained from submillimeter con- 
tinuum data and vary strongly with radius. Figure^shows that 
freeze-out can - for a given age - only occur for a restricted 
region of density and temperature in the envelope. At high tem- 
peratures T > Tev the molecule evaporates whereas at low den- 
sities n < «de the freeze-out timescale is too long. Due to the 
exponential dependence in Eq. Q, the thermal desorption pro- 
ceeds very rapidly as soon as the temperature is higher than 
Tev. In contrast, the freeze-out timescale varies more slowly 
with depth in the envelope, due to the inverse dependence on 
density. 

In order to quantify this scenario, a "drop" abundance struc- 
ture is introduced as a trial profile, with a depleted abundance 
Xd where T < T^^ and n > n^e, and an undepleted abundance 
Xo where T > or n < nje- Whereas Tgy is in principle 
a well-defined quantity depending primarily on the ice mantle 
properties, «de depends on the lifetime of the core compared 
to the depletion timescale (induced in the pre- and protostellar 
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stages) and the dynamical evolution of the core at earlier stages, 
including replenishment of undepleted material through infall 
from larger radii. This may move undepleted envelope mate- 
rial to smaller radii, or higher densities, increasing the derived 
value of nde- On the other h and, within the standard inside-out 
collapse model the infall radii for these very early 

stages are so small that the bulk of the envelope is not yet in- 
falling, including the material near the depletion radius where 
n = ride. 

For this discussion we only assume step functions 
at Tev and «de, which are thus free parameters together 
with Xq and Xu- With these assumptions the full Monte 
Carlo line radiative transfer is performed as described in 
[j0r2ensen et al. (2002) and Schoier et al. (2002) using the code 
of iHogerheiide & van der TakI ( 120001) . We adopt the molec- 
ular data summarized in the Leiden Atomic and Molecular 
Database ( Schoier et al., 2005)'. Single-dish da ta on CO and 
HCO^ taken from Jcirgensen et ZI (l2002l E004hi) are fitted ini- 
tially for four objects at different stages of evolution, leaving 
all 4 parameters free in the modeling of the CO data. The abun- 
dance of HCO^ is expected to reflect the freeze-out of CO as 
illustrated by ch emical models (e.g.. iBergin & Lang e 
lLeeetalll2004 . but the HCO+ lines have higher critical den- 
sities and therefore provide an independent probe. iLee et alJ 
modeled the chemical and physical evolution of a collapsing 
core through its pre- and protostellar stages and found simi- 
lar "drop profiles" for the CO abundances and closely related 
species. For the HCO^ fits, Tgy and Wde are therefore taken from 
the fits to the CO data, giving low values of the ;t'^-estimator. 
With Hde constrained from fits to the observed CO lines, Eq. ^ 
then gives the depletion timescale t^e - IjA for a source with 
a given temperature and density profile. Subsequently, the CO 
data for the entire sample of 16 sources have been fitted, keep- 
ing only Xd and nde as free parameters. 

The pre-stellar cores require a special discussion because 
they do not have a central source of heating so that there is 
only a single freeze-out radius /?de corresponding to density nde- 
The precise value of nde depends in this case on the adopted 
physical structure and abundance profile. As an example of 
the uncertainties, our best-fit models for L1544 are compared 
with other recent work by Bacmann et al. (2002), Tafalla et al. 
( I2OO2I) and iLee et alJ ( 120031) in Fig. El which shows the best 
fit physical structures together with the derived CO abundance 
structures. There are some differences in the adopted under- 
lying density profi les and derived parameters. For example, 
iTafalla et al.l (I2OO2I) use an exponentially decreasing CO abun- 
dance structure toward the center X(r) - Xq exp(-n(r)/nQ) with 
no - 5.5 X lO"* cm"^. The origin of the difference between 
their no and the lower nde of 1 .5 x 10"^ cm'-^ derived using 
step functions of the CO abundance dLee et all l2003[ this pa- 
pe r) stems from (i) th e slightly low er density profile derived 
by lEvans etaP ( 1200 ll) and used by iLee et alJ and in this pa- 
per; (ii) the higher (u ndepleted) CO abundance in the outer re- 
gion of the core from lLee et alJ and this paper, and ( iii) the fa ct 
that CO disappears complet ely in the m odels o f Tafalla et alJ at 
radii less than ~ 5000 AU. iLee et all test different abundance 

' http : //www. strw. leidenuniv.nl/~moldata 
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Fig. 2. Comparison between t he density (bla c k) an d CO 
abundance profiles (grey) fro m iBacmann e t al. (2002) (dot- 
ted lines), iTafalla et all i2002h (solid lines), .Lee et aL C2003l) 
and this paper (dashed lines) for the pre-stellar core L1544. 
The dashed-dot t ed line indicates the constant abundance from 
|j0rgensen et al.l (l2002h . The freeze-out radius at which a sig- 
nificant d rop in the CO a bundance occurs in the step-function 
model of iLee et'all (^2003) is indicated by the vertic al arrow . 
The o pen and solid stars indicate the values of no (.Tafalla et all 
l2002h and nde (this work), respectively. See text for details. 

profile types and find that a step function with a characteris- 
tic density marked with the vert ical arrow in Fig. |2] provides 
the best fit. lBacmann et alj(l2002l) derive only an average abun- 
dance toward the center of the core with a density profile with a 
less dense inner flattened region, and therefore do not give con- 
straints on nde. Taking this difference in the underlying density 
profiles into account their average (or constant) CO abundance 
of 7x10"^ is in good agreeme nt with the constant abu ndance of 
5x10^^ for L1544 derived by |j0rgensen et alJ ( l2002l) . 

3. Results 

The best fit parameters to the CO and HCO^ data for the 4 

sources are listed in Table [2 and illustrated in Fig. |4j|5l For 

the class I object, TMRl, the data are found to be consistent 

with a model without a depletion zone. For CO, Xq is found 

to be 2.7x10""* consistent with t he direct measure ment of the 

CO abundance relative to H? bv lLacv etal](ll994 m a warm 

cloud in which CO is not fr ozen out. Th is is also the expected 

value based on the results of l.T0r^ensen et al. (2002) who found 

a maximum constant CO abundance of ~2xl0"^. Van der Tak 
• 1 

et al. (2000) found a similar maximum constant abundance for 
a sample of high-mass YSOs. In the context of the drop abun- 
dance model presented in this paper, this constant abundance 
applies to sources such as TMRl, where the drop region is so 
small that the abundance profile can be assumed to be essen- 
tially constant. 

For the class objects, in contrast, the constant abundances 
found by J0rgensen et al- (200Z) are about an order of magni- 
tude lower, indicating that the region of CO depletion (the drop 
zone) is large for these objects. These are also the objects where 
the intensities of the low excitation 7=1-0 transitions are un- 
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Table 1. Best fit parameters for CO and HCO^ lines for the pro- 

tostars N1333-I2, L723, TMRl and the pre-stellar core L1544. c^'o (i-o; oso) ci=o (e-i; jcmt) c^'o (3-2; jcmt) 
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N 1333-12 


L723 


TMRl 


LI 544 


Tev [K] 


> 35 


40 






Wde [10* cm-3] 


7 


4 




1.5 












^0 [lO-'*] 


2.7 


2.7 


2.7 


2.7 


[10-"] 


<0.14 


0.14 




0.02 


X'/N 


4.4/7 


2.4/6 


0.75 / 6 


4.2/4 


HCO+: 










Xo [lO-**] 


1.8 


<2.0 


2.7* 




Xu no-'] 


0.26 


0.35 






X^IN 


4.1/5 


1.8/3 


2.1*/ 2 




?de [10^ yrs]" 


2 


3 


< 1" 


8 



Notes: The abundances have been derived from observations of the 
optically thin isotopic species, C' ^O and H'^CO^ assum - 

ing the standard isotopic ratios adopted in I.Tc)rgensen et alJ i2004lll . 
"Depletion timescale corresponding to the derived Wde using Eq. J2j. 
''Model including excitation through collisions with electrons at den- 
sities < 3 X lO'* cm"^. "^Icr limit assuming T^^ of 35 K. 
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Fig. 3. Fitted C^^O line-profiles for NGC 1333-IRAS2. Upper 
panels: constant fractional abundance [^^CO] of 2.4x10"^ from 
|j0r2ensenet al. (2002). Lower panel: best fit drop abundance 
model with Xq = 2.1 x 10 and = 1.4 x 10"^ from this pa- 
per The calibrational uncertainty for each observed spectrum 
is about 20%. 



derestimated unless the drop abundance structure is introduced. 
This is clearly illustrated in Fig. |3] which shows the best fit to 
the C'^O lines fo r NGC 1333-1RAS2 w ith a constant abun- 
dance model from lj0rgensen etaD ( l2002l) and the drop model 
from this paper A si milar plot ill ustrating the fits for L723 can 
be found in Fig. 7 of |j0rgensen et al. ( 2004b) . Fig.l^illustrates 
the confidence levels for the various parameters for NGC 1333- 
IRAS2. It is seen that Tev only has a lower limit and Xd only an 
upper Umit: for Tev higher than 35-40 K the innermost region 
with high CO abundances has a very low filling factor in the 
single-dish beam and therefore only contributes a small frac- 
tion to the beam-averaged column densities. The best fit value 
of X{) is again 2.5-3x10"'*. 

For all sources listed in Table^except TMRl, is found 
to be typically an order of magnitude lower than Xq for both CO 
and HCO^. The CO evaporation temperatures r^y are all found 
to be > 35 K, consistent with the results of lj0rgensen et alj 
(l2002l) . whereas the values of are all of order 2-3x10 with 
a 10-20% statistical uncertainty. 

The derived CO abundance profile depends somewhat on 
the adopted outer radius and any contribution from the sur- 
rounding cloud to the lowest excitation lines. However, the 
same drop abundance profiles are also necessary to reproduce 
higher angular resolution interferometer data, where contribu- 
tions from the larger scale cloud are resolved out, as well as 
to model transitions with higher critical densities of H2CO 
and HCO^ which are not excited in the suiTounding cloud 
JSchoier et alll2OO4lj0rgensenll2OO4 and TableQ. 

Given the success of these models, the C'^0 and C'^O 
lines for the entire sample of J0rgensen et al- (.2002) have been 
fitted to infer the CO abundance structure and the depletion 
time scales. Based on the results for the above 4 sources, it 
seems reasonable to take the undepleted abundance Xq fixed at 



2.7x 10 and the evaporation temperature Tev at 35 K to reduce 
the number of free parameters. We stress that uncertainties in 
the underlying physical mod el such as due to the unc ertain dust 
properties (see discussion in lj0rgensen et al.L l2002h and other 
factors may result in systematic uncertainties by factors of 2- 
3 in the derived absolute abundances and values of n^e and t^^ 
(see also above example on L1544). However, the conclusion 
that a drop abundance profile provides a much better fit than a 
constant abundance profile is not affected by these uncertain- 
ties. Also, since all sources are analyzed in the same way, the 
relative values and inferred trends should be more reliable. 

The fitted Wde and Xd, together with the value and 
derived fje are given in Table |5| For each source, all lines 
could be well fitted ix^^^i ^ 2) using depletion densities 
Hjjg =lxl0''-6xl0^ cm"-' (lO^^""* cm"-') corresponding to de- 
pletion timescales of 2 x lO'^-S x 10^ years (10^=^°-^ years). No 
significant difference in timescales between Class and 1 ob- 
jects is found. For the objects with the least massive envelopes 
{M <, OA Mq) where photodesorption could become an issue, 
the derived fde is a lower limit to the actual depletion timescale. 
In general, however, the depletion radius is located further in 
the envelope than the radius where Ay = 2 for our objects. 

4. Discussion 

The fact that the drop abundance profiles provide better fits than 
constant abundance profiles for all sources and molecules stud- 
ied shows that this structure provides a good representation of 
the dominant chemistry in protostellar envelopes, with other 
chemical effects being of secondary importance. To illustrate 
the dependence on envelope mass. Fig. |6l compares the radii 
corresponding to T^^ and «de for varying envelope masses for 
a (typical) 3 central source. For an object with a large en- 
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Fig. 4. confidence plots for fits to the CO lines toward 
NGC 1333-IRAS2. In each plot the grey area marks the Icr 
confidence region and the subsequent Line contours the 2cr, 3cr 
and 4cr confidence regions. 




Fig. 5. As Fig.lHbut for TMRl (left) and L1544 (right). In the 
left panel the dashed region indicates constant abundance mod- 
els (i.e., models where the radius corresponding to «de is lo- 
cated closer to the central source than the radius corresponding 

to Tev). 



velope mass like NGC 1333-IRAS2 (1.7 Mq), the depletion 
zone is a few thousand AU, a substantial fraction of the en- 
tire envelope. For an object with a small envelope mass like 
TMRl (0.1 Mq), the freeze-out radius moves inwards and the 
size of the depleted region becomes vanishingly small. This in- 
dicates that the trend of increasing (constant) abundances with 
decreasing e nvelope mass for these species ( J0rge nsen et al-l 
reflects the size of the depletion region. 

Based on the above fits, we propose a chemical structure 
of the outer envelopes as shown schematically in Fig. The 
main difference between the pre- and protostellar cores is the 
inner source of heating in the protostars that causes CO to be 
evaporated rather than depleted toward the source center. The 
main difference between the Class and I objects is the size of 
the depletion zone. 

Is this chemical structure also an evolutionary indicator? 
Due to progressive dispersion of the envelope by outflows 
and mass accretion, an intially massive envelope like that of 
NGC 1333-IRAS2 should eventually go through a phase with 



Table 2. Values of Wd and constrained from the CO lines for 
each individual source in the sample, together with the total 
number of fitted lines and derived fde. 



Source 


[cm"'] 






y'- 1 N 

A / ^ ' 


fd 

fyrl 


L1448-I2 


4.8x10^ 


3.0x10" 


-7 


4.2/4 


3x10" 


L1448-C 


6.0x10'' 


2.0x10" 


-5 


10.7 /7 


2x10' 


N1333-I2 


7.0x10" 


1.4x10" 


-5 


4.4/7 


2x10' 


N1333-I4A 


6.0x10^ 


2.0x10" 


-7 


11.2/7 


2x10" 


N1333-I4B 


1.7x10' 


4.3x10" 


-6 


9.5 /7 


7x10" 


L1527 


1.4x10' 


2.7x10" 


-5 


9211 


1x10' 


VT A 1 623" 












L483'' 


1.5x10' 


5.0x10" 


-6 


1.5/3 


9x10" 


L723 


4.0x10" 


1.4x10" 


-5 


2.4/6 


3x10' 


L1157 


<2.3xlO' 








>5xlO" 


CB244 


1.3x10' 


1.6x10" 


-5 


3.8/5 


1x10' 


L1489 


7.0x10'' 


1.0x10" 


-5 


8.1/6 


1x10' 


TMRl 


> 1.0x10' 






0.75 / 6 


<8xlO" 


LI 544 


1.5x10'' 


1.6x10" 


-6 


4.2/4 


8x10' 


L1689B 


> 1.0x10" 


<5xlO- 


5 


0.16/3 


<lxlO** 


IRAS 16293" 


1.0x10' 








1x10' 



"For VLA1623 the 1-0 line observations are significantly underes- 
timated even for a high constant abundance. This likely reflects the 
dense ridge of material in which this source is locate d, which also af- 
fects the higher excitation lines (see also discussion in |j0rgensen et all 
l2002h . '' From fits to single-dish and interferometric C'^O obser- 
vations jj ^^nsej. [2004). 'From fits to H2CO interferometer data 
iSchoier et al.', ■20041.'The CO isotopic lines for IRAS 16293-2422 
are consistent with a constant abundance throughout the envelope but 
do not rule out a drop in abundance. 



< 



10000 



1000 



100 



Lum = 3 Lq 


t - 10^ years - 

de ^ 


. depletion 
zone 


.-■>''''' t - 10'*^ years 








~ " T" "=" 35 K 

ev 



0.1 



1.0 

Mass / Mq 



10.0 



Fig. 6. Characteristic radii as functions of envelope mass for 
a 3 Lq protostar using the models of J0rsensenet al. (2 0o3) . 
An n oc r '-^ envelope density profile was assumed with inner 
and outer radii of 25 and 15000 AU, respectively. The dashed 
line indicate the radius where Tev = 35 K and the solid and 
dotted lines the radii where the depletion timescales are 10'* 
and 10^ years, respectively. The dashed-dotted line indicates 
the difference between the depletion and evaporation radii (i.e., 
size of the depletion zone) for a depletion timescale of 10'' 
years. 
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Fig. 7. Proposed chemical structure for low-mass pre- and pro- 
tostellar objects. The left column gives the temperature and 
density as functions of radius (black solid and grey dashed 
lines, respectively) for three archetypical low-mass pre- and 
protostellar objects: L1544 (pre-stellar core), N1333-I2 (class 
0, Menv > 0.5 Mq protostai-) and TMRl (class I, Mgnv < 0.5 Mq 
protostar). The black dotted lines indicate the derived abun- 
dance structures. The right column gives the depletion signa- 
ture for each class of object with, going from the outside to the 
inside, the dark grey indicating the region where the density is 
too low for depletion (n < nje), the black indicating the region 
where the molecules deplete and the light grey indicating the 
region where they evaporate (T > T^y). 

a low-mass envelope like that around TMRl. However, in con- 
trast with the size of the depletion zone, no clear correlation is 
seen between the derived age f(je and derived envelope mass or 
luminosity for the entire sample (e.g.. Table 1). This indicates 
that the depletion timescales must also reflect other properties, 
such as the mass of the core from which the protostar is formed. 
In other words, an other object like TMRl did not necessarily 
start out with an envelope as massive as that of NGC 1333- 
IRAS2. The relatively short timescales of 10^^*'^ years sug- 
gest that the depletion structure is established in the pre-stellar 
stages, but only after the pre-stellar core has become dense 
enough that freeze-out really sets in. The rather small scatter 
seen indicates that this stage is short, of order 10^ years. 

This timescale provides an interesting independent age con- 
straint for the studies of core collapse, since it is shorter than 
the a ges derived from the above mentioned statistical studie s 
(e.g.. iLee & MversL Il999t Ijessop & Ward -Thompson' '2000"). 
This apparent discrepancy may reflect simply the definition 
of the pre-stellar stage: our estimate refers only to the dense 
pre-stellar phase where depletion has become significant. On 
the other hand, the statistical studies may be missing low- 



luminosity embedded infr ared sources due to the limited sen- 
sitivty of IRAS (see, e.g. jYoung et al.Ll2004l) and thus lead to 
overestimates of the timescales. Further unbiased surveys, e.g., 
with the Spitzer Space Telescope, may shed further light into 
this issue. 

To refine the proposed models, it will be necessary to cou- 
ple mod els for the dynamical and radial chemical evolution as 
done bv lLeeet a\\ h()()4) to fully address the use of «de as a 
tracer of age. It is interesting to note that our derived empirical 
drop abundance structure agrees well with these detailed mod- 
els, not only qualitatively but also quantitatively. As their Fig. 6 
shows, the timescales over which the heavy depletion occurs in 
the pre- and protostellar stages is only of order 10^ years. 

Another important consequence of the "drop" chemical 
structure is that it affects the tracers of infall in protostostellar 
envelopes. Since the dynamical structure of protostellar cores 
is often infe rred from fits of line pro files of molecules such as 
HCO^ (e.g.. iGregersen et al.L Il997l) . knowledge about the ra- 
dial chemical structure is important for detailed descriptions of 
the infalling envelope. For example, the location of the collapse 
radius in the inside- out collapse model for th e envelope around 
NGC 1333-IRAS2 il^rgensen etani2004ah is at =^ 1000 AU 
which is in the middle of the drop zone where the temperature 
is w 25 K. The exact "infall" line profile will therefore depend 
critically not only on the velocity field but also the presence 
and location of the outer pre-depletion {n < «de) zone and the 
amount of depletion. 
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